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In present talk I show that primordial helical turbulence produced during a first-order 
phase transition induces circularly polarized cosmological gravitational waves (GWs). 
The degree of polarization as well as the characteristic frequency and the amplitude 
of these GWs depend crucially on phase-transition model of primordial turbulence. I 
present a brief discussion on the possibility of detection. 

1. Introduction 

There are several astrophysical observations indicating that the magnetic field in 
the Sun, and some galaxies and clusters of galaxies might have an helical structure 
^. On the other hand, recently it has been realized that magnetic (or/and kinetic) 
helicity can be generated during an early epoch of the universe Primordial 
helical motions influence magneto-hydrodynamical processes in the early plasma 
as well as cosmological perturbation dynamics Since cosmological gravitational 
waves (GWs) propagate without significant interaction after they are produced, 
once detected they should provide a powerful tool for studying the early Universe 
at the time of GW generation ^. For the case at hand, GWs are considered as a 
test for primordial helical turbulence. The present formalism is general and can be 
applied to study the generation of stochastic GWs by any helical vector field (e.g., 
helical magnetic fields "^'S'^). 

It should be noted that various mechanisms for GW generation have been stud- 
ied, including: quantum fluctuations during inflation bubble wall motion and 
collisions during phase transitions cosmological magnetic flelds i2,i3. g^j^^^ plasma 
turbulence I3'i4,i5^ Here I focus on gravitational radiation generated by primordial 
helical turbulence This talk is based on results obtained in collaboration with G. 
Gogoberidze, A. Kosowsky, A. Mack, and B. Ratra (see Refs. i3,i6-j^ g^^^ thaX he- 
lical turbulence generates circularly polarized stochastic gravitational waves (GWs) 
and we compute the polarization degree. Primordial polarized GWs might be gen- 
erated from quantum fluctuations accounting for the gravitational Chern-Simons 
term ^'^ . 

As it is well known GWs are generated by the transverse and traceless part 
of the stress-energy tensor T^^ . In our case describes a turbulent cos- 
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mological fluid after a phase transition iiii3,i4,i5_ pgj. spatial indices i j, 
Ty (x) = {p + p)ui{x)uj{x), where p and p are the fluid pressure and energy density 
and u(x) is the fluid velocity. The fluid enthalpy density p + p is taken to be con- 
stant throughout space. The transverse and traceless part of T^j in Fourier space is 
Ilij{k,t) = [Pii{k)Pjm{'k) - iPij(k)Pi„(k)]r;„(k,f), where the projector onto the 
transverse plane Pij(k) = Sij — kikj with ki = ki/k, where ki is the wave-vector. 
Consistent with observations we have assumed flat space sections. 

Since turbulence has a stochastic nature the induced GWs should be also 
stochastic. They will form relic GWs background, and if there is a parity breaking 
in the source, one can expect the circular polarization of such waves To compute 
the gravitational waves signal from helical turbulence, first the model of turbulence 
should be specified. 

2. Model helical turbulence 

To model the turbulence we assume that in the early Universe at time tin (at a 
phase transition) liberated vacuum energy pvac is converted into (turbulent) kinetic 

energy of the cosmological plasma with an efficiency k over a time scale TsUi on 
a characteristic source length scale Ls After generation, the turbulence kinetic 
energy cascades from larger to smaller scales. The cascade stops at a damping scale, 
Lu, where the turbulence energy is removed by dissipation. As usual, we assume that 
the turbulence is produced in a time much less than the Hubble time, Tgtir ^ 1 /-ffin 
— here iJin is the Hubble parameter at t-m and rapidly generates GWs. We 

therefore ignore the expansion of the Universe when studying the generation of GWs 
(see below). 

To compute the induced GW power spectrum one must have the source two- 
point function (H* (k, i)H;„i(k', t')). This is determined by the fluid velocity two- 
point correlation function. For stationary, isotropic and homogeneous flow the ve- 
locity two-point function is ^'^ 

«(k)w,(k')) = {2nf6^^\k - k')[PijPs{k) + icijikPHik)]. (1) 

Here Ps{k) and Pnik) are the symmetric (related to the kinetic energy density 
per unit enthalpy of the fluid) and helical (related to the average kinetic helicity 
(u • (V X u))) parts of the velocity power spectrum Causality requires Psik) > 
\PH{k)\ (Schwartz inequality, see p. 161 of Ref. l^ and Refs. 20,21,2,3,4) 

However, in the present case, the source of turbulence acts for only a short 
time T, possibly not exceeding the large-scale eddy turnover time rg (corresponding 
to length scale Ls) — for self-consistency, however, we assume that the source is 
active over a time t — max^Tsth , ts) resulting in a time-dependent velocity 

spectrum. To model the development of helical turbulence during the time interval 
{tin,tfi = tin + t) we make several simplifying assumptions: 

(a) Turbulent fluid kinetic energy is present on all scales in the inertial (Kol- 
mogorov) range ks < k < ko- Here ks = ^tt/Ls and kn = I-k/Ld, and the inertial 
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range includes length scales L G [Ld,Ls), i.e., length scales shorter than those 
on which energy flows into the turbulence and larger than those on which the tur- 
bulence energy is dissipated. We also assume that the energy is injected into the 
turbulence continuously over a time r, rather than as an instantaneous impulse 

13,14,15,16 

(b) Unequal time correlations are modeled as ^^'^^ 

«(k,i)u,(k',t')> = {2T:f5^^\\ii-'k')[PijFs{k,t-t')+ieijikiFH{k,t-t% (2) 

where the t — t' dependence of the functions Fs and Fh reflects the assumption 
of time translation invariance. Since energy is injected continuously, at t = t' € 
(<in,te), Fs{k,0) = Ps(k) and FH(k,0) = Pnik). 

(c) The decay of non-helical turbulence is determined by a monotonically de- 
creasing function Di{t) and Fs{k,t) = Ps{k)Di{t), p. 259 of Ref. Extending this 
assumption to the helical turbulence case we also model FH{k,t) = PH{k)D2{t), 
where D2{t) is another monotonically decreasing function. Since in the considered 
model most of the power is in the inertial range for simplicity we discard power 
outside the inertial range by truncating Ps and Ph at k < ks and k > kn, i.e., we 
assume Ps and Ph (and so Fs and Fh) vanish outside the inertial range. 

(d) We model the power spectra by power laws, Ps{k) oc A:"^ and Pnik) oc 
fc"". For non- helical hydro dynamical turbulence the Kolmogorov spectrum has 
ns = —11/3. It has been speculated that in a magnetized medium an Iroshnikov- 
Kraichnan spectrum with ns = —7/2 might develop instead. The presence of hydro- 
dynamical helicity makes the situation more complex. Two possibilities have been 
discussed. First, with a forward cascade (from large to small scales) of both energy 
and helicity (dominated by energy dissipation on small scales) one has spectral in- 
dices ns = —11/3 and nn = —14/3 (the helical Kolmogorov (HK) spectrum), p. 243 
of Ref. Second, if helicity transfer and small-scale helicity dissipation dominate, 
ns = 7iH = —13/3 (the helicity transfer (HT) spectrum) The HK spectrum 
has been observed in the inertial range of weakly helical turbulence (i.e., where 
\PH{k)\ Ps{k)) For strongly helical hydrodynamical turbulence the charac- 
teristic small-scale length scale of helicity dissipation is larger than the Kolmogorov 
energy dissipation length scale Therefore here the inertial range is taken to con- 
sist of two sub-intervals, both with power-law spectra. For smaller k the spectra 
are determined by helicity transfer and have ns = tih = —13/3, while for larger 
k turbulence becomes non-helical and the more common HK spectrum is realized. 
Since GW generation is mostly determined by the physics at small k (near ks) 

it is fair to only use the HT spectrum in this case also 

Based on these considerations we model the primordial spectra as Ps{k) = 
S'ofc"^ and Pnik) = AoA;^^"""^"", where: (i) for the HK case So = vr^Cfee^/^ 
and Aq = 7r^Cfe5/(£^/^A:s) implying ^o/5'o = S/{eks); and, (ii) for the HT case 

= -K^CsS^I^ and = 'n^CaP'l'^ Here e and Z are the energy and mean helicity 
dissipation rates per unit enthalpy, and Cfc (the Kolmogorov constant), Cg, and Ca 
are constants of order unity. 
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The formalism described here is appUcable for locahy isotropic turbulence. An- 
other requirement is connected with plasma viscosity, i.e. the Reynolds number Re, 
which is related to the cut-off scales ks and ko (for HK case, Rc = (ko/ksy^^), 
has to be enough large '^^ . Summarizing, to model cosmological turbulence it is 
required to know (i) the part of the vacuum energy converted to the turbulent 
motions KPvac; (ii) the "stirring" scale Ls; (iii) the temperature of the universe Tin 
when turbulence is generated. This temperature defines the energy density and en- 
thalpy {p + p) of the plasma; (iv) the ratio between energy and helicity dissipation 
rates S/{kse). Knowing these quantities and the turbulent motion spectral index, 
we are able to obtain the damping scale ko and the plasma viscosity v (p. 483 of 
Rcf. ^^). In particular, for locally isotropic turbulence the energy dissipation rate 
£ = 21/ Ji^^ dk k^Ps{k)/Tr'^ is equal to the source power input, i.e., 

U=^, (3) 

3 pT 

where k is the phase transition efficiency. For HK turbulence with ng = —11/3, the 
damping scale is related to plasma viscosity via 

Ij^^rkj, = (4) 

For the HT spectrum {ns = —13/3) the damping scale is connected with plasma 
viscosity and the dissipation rates e and S via 3i'(5kD)^^^ = £; using the fact that 
the total energy density of turbulent motion pturb = ^Pvao we find 



K 



ks_ 
ko 



2/3 



2,Tvki (5) 



The L-scale eddy tiirnover (circulation) time can be obtained as a ratio of length 
scale L to the physical velocity vl related to the normalization of the symmetric 
part of turbulent motion spectrum Ps{k) (or the energy dissipation rate s). It is 
easy to find that 2tl — 2>e~^l^l?'l^ for any L inside the inertial range Ljj < L < Lg 

13,16 

It should be underlined that we use the non-rclativistic turbulence model to 
describe processes in the early universe. The Kolmogorov model of turbulence has 
been proved only for non-relativistic velocities in plasma, and there is no standard 
model for the relativistic case. We have also used the non-compressible plasma 
approximation. For the realistic case of a compressible fiuid one can expect that a 
large enough input of energy may lead to the formation of shocks. On the other hand, 
our assumption that the upper limit of plasma velocity is equal to sound speed may 
be justified because of significant thermal dissipation due to shock waves. Presuming 
that shock fronts accumulate kinetic energy, the gravitational radiation signal will 
be increased (relative to the non-relativistic case), if we account for relativistic 
effects 13.14. 
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3. gravitational wave polarization degree 

Since the turbulent source acts a short time (comparing with Hubble time-scale), 
the expansion of the universe can be neglected in the GW equation of motion, then 
in wave-number space we get 



hij (k, t) + k'^hij (k, t) = le-irGUij (k, t). (6) 

Here G is the Newtonian gravitational constant. We use natural units h = 1 = c, 
physical/proper wave- numbers (not co- moving ones), and an over-dot denotes a 
derivative with respect to time t. /ijj(k, t) = J (fxe^^'^hij{x,t) and hij{x,t) = 
J cfik e~^^'^ hij (k, t) / {2tt)^ is the Fourier transform pair of the transvcrse-traceless 
tensor metric perturbation which is defined in the terms of the complete metric 
perturbation hj j = Sgij subject to the conditions hu = and /ijjfc-' = 0. Choosing 
the coordinate system so that unit vector 63 points in the GW propagation direction, 
using the usual circular polarization basis tensors efj = — (ei±ie2)i x {ei±ie2)j/V^ 

and defining two states h'^ and h~ corresponding to right- and left-handed 
circularly polarized GWs, we have hij = h'^efj + h~e~j. 

Stochastic turbulent fluctuations generate stochastic GWs. Gaussian-distributed 
GWs may be characterized by the wavenumber-space two-point function 

{h*j{k,t)him{i^',t)} = (27r)35(3)(k _ k') [MijlmH{k,t) + iA,7r„H(A;, i)] . (7) 

Here H{k,t) and 'H{k,t) characterize the GW amplitude and polarization, 

4A^ij777i(k) = PilPjm -t- PimPjl PijPlm-} and S^jj/m (k) = \^q(^Pj^€ilq -\- Pil^^jmq "t~ 

Pim^jiq + Pji^imq) are tensors, and Ciji is the 3 dimentional fully antisymmetric 
symbol. 

The GW degree of circular polarization is given by 

^ {h+*{k,t)h+{k',t) - h-*{k,t)h-{k',t)) ^ H(fc,t) 

{h+*{k,t)h+{k',t) + h-*{k,t)h-{k',t)) H{k,ty ^' 

Both H{k,t) and H{k,t) are obtained by solving Eq. (6), and are related to 
Ilij{k,t). For instance, an axisymmetric stochastic vector source (non- helical tur- 
bulent motion or any other non-helical vector field) induces unpolarizcd GWs with 
|/i"'"(k,t)| = |/i~(k,f)| I2,i3,i4_ other hand, the presence of a helical source 

(i.e., a source that breaks parity symmetry) alters this situation 

Given a model of the turbulence, the turbulent source two-point function is 

{U*j{k,t)Uim{k',t + y)) = {2wf6^''\k-k') [Mijimf{k,y) + iAijimg{k,y)] , (9) 
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where Mijim and Aijim are defined below Eq. (7). The functions f{k, y) and g{k, y) 
that describe the symmetric and helical parts of the two-point function are 

/(fc, y) = / d'Pij d'p2S^'^ (k - - [(1 + 7^) (1+ 

+ (3')Dj{y)Ps{pi)Ps{p2) + ^iPDliy) ] Ph{pi)Ph{p2), (10) 
9{k, y) = ^P + Pf^^^)^-^y) J d^p,J ,3,, ,(3) (k _ _ p,) [(1+ 

+ 7')/3Ps(pi)Ph(P2)(1 + P^hPH{pi)Ps{P2) ] , (11) 

where 7 = k - pi and /3 = k - P2. The antisymmetric (parity-odd) source term g{k, y) 
vanishes for turbulence without helicity. 

To determine II{k, t) and W(fc, t) we solve Eq. (6) assuming that there is no GW 

for times t < tin, i-e., we choose as initial conditions ft,ij(k,tin) = = hij{k.,ti^). 
To compute the induced GW power spectrum we use the averaging technique de- 
scribed in detail in Sec. III.B of Ref. The main points are: (i) the (5^^'(k — k') 
in Eqs. (7) and (9) ensure statistical isotropy of the GWs and allow k' and k to 
be interchanged; (ii) the statistical average can be approximated by either a time 
or a space average (this is justified for locally isotropic turbulence where the cor- 
relation function (wj(xo, fo)wj(xo +r,to +t)) for t < r does not depend on to for 
relatively small t, sec Sec. 21.2 of Rcf. ^'''); (iii) we choose to time average since the 
Green function for Eq. (6) and the source term Him (k, t) are time dependent (we 
then have three time integrals for the GW two-point function, two from the GW 
Green function solutions and one from the time averaging, J*^'^'^ dt/T, whore ti 
is an arbitrary time during the time interval r when the source is active and the 
averaging time T < r); (iv) we take Tk:^ 1 (since k^^ is of order the light crossing 
time for scale L k~^ while T is of order the L scale eddy turnover time). These 
approximations result in Eq. (32) of Ref. 

(/i*.(k, ffi)/ii„(k', ifi)) ~ ^^1^7^ dt cos(fci)(n*.(k, ii)ni„(k', h + 1)). (12) 

The two-point function on the r.h.s. of the integral is independent of ti and 
(k, tfi)/i;m(k', ffi)) is proportional to the source duration time r, as expected 
for locally isotropic turbulence (for a more detailed discussion see p. 358 of Ref. ^^). 

Prom Eqs. (9)-(12) we see that the symmetric H{k,t) and helical H{k,t) parts 
of the GW two-point function in Eq. (7) are integrals over y of cos(A;y)Z?f (y), 
cos{ky)D2{y), and cos{ky)Di{y)D2{y). Both Di{y) and £^2(2/) are positive 
monotonically- decreasing functions of y, and since cos(fcy) oscillates on a time scale 
shorter than the characteristic decay time for the D(?/)'s, the integrands oscillate 
and dy cos{ky)Fa{p, j/)Fb(|k - p|, 2/) ~ Pa{p)Pb{\k - p\)/{V2k) (where a and b 
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can be S oi H). Integrating over angles, we find at t = t{i 

H{k) c^AjdpipiJ dp2 P2e [(1 + 7^)(1 + Pl)Ps{pi)Ps{p2)+ 

+ 4jp/3j,Ph{pi)Ph{P2) ] , (13) 

Hik) ^2aJ dpi pi J dp2 P2Q [{1+iI)PpPs{pi)Ph{p2) + 

+ (l^j3l)^^PH{pi)Ps{p2)]. (14) 

Here A = ar/iiTr^k^) where a = V2ip + pfiSnGf , 7p + pj -pl)/{2kpi), 

Pp = {P+pl-pl)/{2kp2), e = 0{pi+P2-k)0{pi + k-p2)0{p2 + k-pi), and 6* is the 
Heaviside step function which is zero (unity) for negative (positive) argument. Let's 
note that present approximations difi'crs from ones used previously in Rcfs. The 
diff'erence comes from using S'^^^ (k— p) function to evaluate angle integrals and don't 
absorb modulus integration |p|. This preserves to make errors when approaching 
the integration edges. In particular , since the spectra -P(|k — p|) are defined only 
for ks < |k — p| < fc/j, then the integration over spatial angle 7 can not range from 
— 1 to 1 as it is usually done making the convolution when the spectra are defined 
in whole of p e (0, 00). 

For power-law Ps{k) oc fc"-^ and Pnik) oc fc"^ the integrals in Eqs. (13) and (14) 
can bo done analytically, but the results arc complicated and do not edify. Instead 
we compute the degree of circular polarization, Eq. (8), by evaluating the integrals 
numerically for different parameter values Results are shown in Fig. 1. 

Figure 1 and other numerical results show that for maximal helicity turbulence 
(when ^0 = 5*0) with equal spectral indices uh = ns < —3, the polarization degree 
V{k) ~ 1 (upper solid line). For weaker helical turbulence (when Aq < So) with 
Uh — ns < —3, V{k) CAq/So, where 1 < C{ns,nH) < 2 is a numerical factor 
that depends on the spectral indices. For HT turbulence with ns = nn = —13/3, 
C ~ 1.50, while for Iroshnikov-Kraichnan MHD turbulence {ns = Uh = —7/2), 
C w 1.39. Excluding the edges of the inertial range kg < k < ko, an analysis of 
Eqs. (13) and (14) shows that the main contribution to the integrals come from 
areas with pi ^ ks, P2 ^ k and pi ~ A;, P2 ~ fcs- In this case (for arbitrary spectral 
indices ns, nn < -3) H{k), n{k) oc fc"s-3fcs^+^ 



4. Results and discussion 

To make the connection with the GW measurements, we have to define the real- 
space two point correlation function (|/i='=(x, fend)P)- We find 

GW amplitudes are conventionally expressed as (/i'-'(x, tfi)/i*^(x, tg)) = 
2/Q°°rfln/[|/i+(/)|2 + \h-{f)W see Eq. (11) of Ref. where the frequency / of 
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Fig. 1. GW polarization degree 'P{K,tfi), Eq. (8), as a function of scaled wave-number K = k/ks 
relative to the large-scale wave-number ks on which energy is pumped into the turbulence. This is 
evaluated at time tg, after the turbulence has switched off, and remains unchanged to the present 
epoch. It has been computed for a damping wave-number ko = lOfcs. Three pairs of curves are 
shown. Solid lines correspond to the amplitude ratio Aq/Sq = 1, dashed lines to Ao/5o = 0.5, and 
dot-daslicd lines are for Aq/So = 0.2. The upper line in each pair corresponds to HT turbulence 
with ns = riH = —13/3 and the lower line to HK turbulence with ng = —11/3 and njj = —14/3. 
Even for helical turbulence with Ao/Sq < 0.5, for large wave-numbers k ~ fcjj, ns = njj = —13/3 
is unlikely so the large K part of the lower dashed and dot-dashed HT curves are unrealistic. The 
large k ^ ko decay of the HK curves is a consequence of vanishing helicity transfer at large k 



a GW generated by an eddy of length L is f = Tj^^ where tl is the eddy turnover 
time 11.13.14 Using / = 2e^/^L-y^/3 we get 



f.' ^. ,s^2L-,^JS^. (16) 

277^ \j \ns + 3\ " y \ns +3| 

where fs is the frequency that corresponds to the stirring length Ls- Both the 
frequency and the amplitude of GWs are inversely proportional to the cosmological 
scale factor, so the frequency / today and / when the temperature was Tin = 
100 Tioo GeV are related by / = 1.65 x lO-^TioosJoo/Z-f^in Hz where c/i„ = 
lOOgioo is the number of relativistic degrees of freedom at ti„. Since we truncate 
turbulence power for L > Ls, the GW spectrum is non-zero only for f > fs, where 

16 



fs = 1.9 X 10-,&±5 C-V^' ( ^V""'^'Vtoo .rt Hr^ Hz (17) 



ns + 31 Vi// \L 



Here we used Eqs.(3)-(5). For HK turbulence {ns = —11/3) the expression for fs 
will be transformed in Eq. (53) of Ref. 

Using Eqs. (15) and (17), and neglecting the weak fc-dependence of the GW 
polarization degree V, we find that K^{f) oc f~^^^^ for the HK case i^.i^^ while for 
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HT turbulence h^{f) oc /^i^^^. We expect such a steeper dependence on frequency 
for helicity induced GWs, since the hehcity transfer rate is more important on 
larger scales. In both cases the amplitude of the GW spectrum peaks at the stirring 
frequency fs- 

To examine the prospect of such circularly polarized GWs detection, we compute 
the energy density parameters of GWs. This depends strongly on the frequency 
band. The GW energy density parameter for frequency /, flowil) is given by (see 
Eq. (7) of Ref ^7) 

^Gw{f)h' = 5.9 X 10'^\h+{f)f + \h-{m{I-f, (18) 

where h is the Hubble constant in units of 100 km sec~^ M-pc~^. In our case, 

o in..^in-ii -i/3 / 4 ng + 5 ^ f Ld\^^"''+''^ f inp^^cLs Y 

naMf)h - 1.05 X 10 gj ■ j j [-^) 

/ I \ 2(2ns+5)/(ns+5) 

As it is expected the stirring frequency fs and the GW spectrum are very sensitive to 
phase transition properties (for details sec Fig. 1 of Ref. ^^). If the phase transition is 
strongly first order, fs — 5x 10^'^Hz (for the HK case) is near the LISA frequency 
range, but the amplitude of the GW signal is below LISA sensitivity i^'^^'^. Also the 
gravitational radiation signal coming from white dwarf binaries will be compatible 
by amplitude with relic signal considered above (for the frequency range around 
10"^ — 10"^ Hz making difficult to distinguish the source of gravitational waves. 
Thus it is unlikely that the circularly-polarized GWs considered here will be detected 
in the near future, however, GWs generated by helical turbulence will have a enough 
high degree of circular polarization and future detector configurations may well 
be able to. 
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